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ABSTRACT 

We present high-speed, three-colour photometry of the echpsing cataclysmic variable 
SDSS J150722.30-t-523039.8 (hereafter SDSS J1507). This system has an orbital period 
of 66.61 minutes, placing it below the observed "period minimum" for cataclysmic 
variables. We determine the system parameters via a parameterised model of the 
I eclipse fitted to the observed lightcurve by minimisation. We obtain a mass ratio 

of 9 = 0.0623 ± 0.0007 and an orbital inclination i = 83°. 63 ± 0°.05. The primary 
mass is — 0.90 ± O.OIA/q. The secondary mass and radius are found to be = 
g ; 0.056 ± O.OOIMq and Rr = 0.096 ± O.OOli?© respectively. We find a distance to the 

• ' system of 160 ± 10 pc. The secondary star in SDSS J1507 has a mass substantially 

OjO I below the hydrogen burning limit, making it the second confirmed sub-stellar donor 

in a cataclysmic variable. 

The very short orbital period of SDSS J1507 is readily explained if the secondary 
star is nuclearly evolved, or if SDSS J1507 formed directly from a detached white 
dwarf/brown dwarf binary. Given the lack of any visible contribution from the sec- 
" ondary star, the very low secondary mass and the low He I A6678/Hq; emission line 

, ratio, we argue that SDSS J1507 probably formed directly from a detached white 

■ dwarf/brown dwarf binary. If confirmed, SDSS J1507 will be the first such system 

identified. The implications for binary star evolution, the brown-dwarf desert and the 
common envelope phase are discussed. 

Key words: binaries: close - binaries: eclipsing - stars: dwarf novae - stars: individual: 
SDSSJ1507-I-5230 - novae, cataclysmic variables 
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1 INTRODUCTION 



Cataclysmic variable stars (CVs) are a class of interacting 
binary system undergoing mass transfer, usually via a gas 
stream and accretion disc, from a Roche-lobe filling sec- 
ondary to a white dwarf primary. A bright spot is formed 
at the intersection of the disc and gas stream, giving rise to 
an 'orbital hump' in the lightcurve at p hases 0.6 — 1.0 due 
to foreshortening of the bright-spot. [Warner! l| 19951 ) gives a 
comprehensive review of CVs. The lightcurves of eclipsing 
CVs can be quite complex, with the accretion disc, white 
dwarf and bright-spot all being eclipsed in rapid succession. 
With sufficient time-resolution, however, this eclipse struc- 



ture allows the system parameters to be determined to a 
high degree of precision (|Wood et alj|l98^ . 

SDSS J1507 was discovered eis part of the systematic ef- 
fort to catalog and characterise the CV s found by the Sloan 
Digital Sky Survey (jSzkodv et aL|[2005l ). Its spectrum shows 
a strong blue continuum and absorption lines from the white 
dwarf, as well as the broad, double-peaked emission lines 
from the accreti on disc, which are characteristic of high incli- 
nation systems llHorne fc Marsh 1986 1. Follow up photome- 
try of SDSS J1507 (|Szkodv et al.ll2005h showed deep eclipses, 
and an orbital period of ~ 67 minutes. The orbital period 
of SDSS J1507 places it below the period minimum, a sharp 
cut-off in the number of CVs at ~ 78 mins. Whilst there are 



2 



Table 2. Mid-eclipse timings (HJD + 2453798). Tlie uncertain- 
ties for the values measured from our eclipses are 8.0xl0~® for 
the red lightcurves and 4.0 xlO"'' for the green lightcurves. 
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objects related to CVs with orbital periods of less than 78 
minutes, known as AM CVn stars, they are believed to con- 
tain a white dwarf and a degenerate heliu m white dwarf or 
a sem i-degenerate helium secondary star l|Nelemans et al.l 
I2OOII ). Only two other CVs exist with orbital periods be - 
low the peri od minimum, V485 Cen lAugusteiin et al.lll996l ) 
and EI Psc (jThorstensen et aLllioO^ ). Possible explanations 
for the presence of CVs below the period minimum include 
a donor which has undergone significant nuclear e volution 
in its core prior to mass transfer (e.g. ,Thorstcn sen et al.l 
l2002h . or the direct formation of a CV from a detached white 
dwarf/brown dwarf binary (|Politanoll2004l ) . The presence of 
deep eclipses in SDSS J1507 allows a precise determination 
of the donor star's mass, which may in turn allow us to 
decide between these two scenarios. 

In this paper we present ultracam u' g'r' lightcurves of 
SDSS J1507, and use these lightcurves to derive the system 
parameters. The observations are described in section [2l the 
results are presented in section [S] and discussed in section [l] 



2 OBSERVATIONS 

On nights between Mar 03''' 2006 and Mar OS*** 2006, 
SDSS J1507 was observed simult aneously in the SDS S-M^ffV 
colour bands using ultracam jDhillon et al.ll2007h on the 
4.2-m William Herschel Telescope (WHT) on La Palma. A 
complete journal of observations is shown in table [1] The 
observations were taken between airmasses of 1.0-1.9, in 
typical seeing conditions of 1.5 arcsecs, but with a range 
of 0.6-5.0 arcsecs. The data were taken in photometric con- 
ditions. Eight eclipses were observed in total. Data reduc- 
tion was carried out in a standard manner usin g the ULTRA- 
CAM p ipeline reduction software, as described in Feline et al.l 
l|2004 ). and a nearby comparison star was used to correct the 



data for transparency variations. Because of the absence of a 
comparison star which was sufficiently bright in the w'-band, 
the w'-band data was corrected using the g'-band data for 
the comparison star, with appropriate corrections for extinc- 
tion. 



3 RESULTS 

3.1 Orbital Ephemeris 

The times of white dwarf mid-ingress Twi and mid-egress 
Two were determined by locating the minimum and maxi- 
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Figure 1. The phase folded u'g'r' (from top to bottom) 
lightcurves of SDSS J1507. All eight observed eclipses are shown 
superimposed to emphasise the lack of variability between them. 



mum times, respectively, of the lightcurve derivative. Mid- 
eclipse times, Tmid, were determined by assuming the white 
dwarf eclipse to be symmetric around phase zero and taking 
Tmid = (Twe + T'wi)/2. Mid-eclipse timings are shown in ta- 
ble [2] The errors on our mid-eclipse times were adjusted to 
give — 1. The g'-band and the r'-band lightcurves only 
were used, given the significantly poorer quality of the u'- 
band lightcurve. The results from the two lightcurves were 
combined with a weighted mean, giving an ephemeris of: 



HJD = 2453798.738844 + 0.04625829 E. 

2 ± 3 

There was no significant deiviation from linearity in the 
O — C times. This ephemeris was used to phase our data for 
the analysis which follows. 



3.2 Lightcurve morphology and variations 

Figure [T] shows the eight observed eclipses of SDSS J1507, 
folded on the orbital phase. The white dwarf ingress and 
egress features are clearly visible, as are the ing ress and 
egres s features of the bright spot (see chapter 2.6.2 of lWarnerl 
1 19951 . for an illustrated example of an eclipse in a typical 
dwarf nova system) . The lightcurve is dominated by emission 
from the bright spot and the white dwarf. The eclipse of the 
accretion disc is difficult to discern by eye, implying that the 
accretion disc contributes little light to this system. There 
is no indication of ellipsoidal variations from the secondary 
star, but a double-humped morphology is seen. This double- 
humped morphology likely arises from the bright spot be- 
ing visible throughout the orbit, which implies the accretion 
disc in this system is optically thin. There is very little vari- 
ation between the individual eclipse lightcurves, implying 
that photometric precision and flickering dominate over sys- 
tematic errors in lightcurve measurement for this system. 
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Table 1. Journal of observations. All observing nights were clear of cloud, but the night of 03'"'' Mar 2006 suffered from high humidity, 
whilst the seeing was poor for the night of OS*** Mar 2006 . The dead-time between exposures was 0.027 s for all observations. The 
relative GPS timestamping on each data point is accurate to 50 fis. 



Date 


Start Phase 


End Phase 


Filters 


Exposure time (s) 


Data points 


Eclipses 


Seeing (arcsec) 


Airmass 


2006 Mar 03 


-0.50 


0.22 


u'gV 


1.995 


1450 


1 


1.1-1.5 


1.092-1.096 


2006 Mar 04 


17.57 


18.18 


u'gV 


1.995 


1206 


1 


1.4-2.8 


1.441-1.630 


2006 Mar 04 


19.83 


20.14 


u'gV 


1.995 


636 


1 


1.4-2.2 


1.138-1.164 


2006 Mar 05 


38.53 


39.16 


u'gV 


1.995 


1257 


1 


2.0-4.0 


1.620-1.900 


2006 Mar 05 


41.39 


42.14 


u'gV 


1.995 


1504 


1 


2.0-5.0 


1.112-1.164 


2006 Mar 07 


86.48 


87.42 


«'gV 


1.995 


1893 


1 


0.8-1.6 


1.096-1.141 


2006 Mar 08 


103.72 


104.29 


u'gV 


1.995 


1161 


1 


1.4-2.4 


1.473-1.666 


2006 Mar 08 


104.91 


105.14 


u'gV 


1.995 


454 


1 


1.0-1.5 


1.286-1.327 



3.3 A parameterised model of the eclipse 

To determine the system parameters we used a physical 
model of the binary system to calculat e eclipse liKh t curves 
for each of the various components. iFeline et al.l (|2004| ) 
showed that this method gives a more robust determination 
of the system parameters in the presence of flickering than 
the derivative method of Wood et al.l ll 19861) . We used the 
technique developed by iHorne et al. (|l994l ) and described 
in detail therein. This model relies on two critical assump- 
tions: that the bright spot lies on the ballistic trajectory 
from the secondary star, and that the white dwarf is accu- 
rately described by a theoretical mass-radius relation. The 
error introduced by the latter assumption can be estimated 
by comparing the results from different theoretical models, 
and we show later in this section that it is smaller than the 
statistical erro rs. Wh i lst th e former assumption cannot be 
directly tested, iFeUnd (|2005l ) show that masses derived with 
this model are consistent with other methods commonly em- 
ployed in cataclysmic variables over a wide range of orbital 
periods, including systems with similar accretion discs to 
SDSS J1507. 

We fit the model to all the observed eclipses, which were 
phase-folded and binned into groups of three data points. 
The data and their model fits are shown in figure [S] The 10 
parameters that control the shape of the lightcurve are as 
follows: 

(i) The mass ratio, q = Mr/M^. 

(ii) The white dwarf eclipse phase full-width at half- 
depth, A(/>. 

(iii) The outer disc radius, Rd/a. 

(iv) The white dwarf limb-darkening coefficient, Uv,. 

(v) The white dwarf radius, Rv,/a. 

(vi) The bright-spot scale, S/a. The bright-spot is mod- 
elled as a linear strip passing through the intersection of 
the gas stream and disc. The intensity distribution along 
this strip is given by {X/S)^e~'^^^ , where X is the distance 
along the strip. 

(vii) The bright-spot tilt angle, Ob, measured relative to 
the line joining the white dwarf and the secondary star. This 
allows adjustment of the phase of the orbital hump. 

(viii) The fraction of bright-spot Ught which is isotropic, 

fiso- 

(ix) The disc exponent, b, describing the power law of the 
radial intensity distribution of the disc. 

(x) A phase offset, (j>o- 



The AMOEBA algorithm (downhill simplex; iPress et al.l 

Il986h was used to adjust all parameters bar Uw to find 
the best fit. A linear regression was used to scale the four 
lightcurves (for the white dwarf, bright-spot, accretion disc 
and secondary) to fit the observed lightcurves in each pass- 
band. The excellent agreement between model and data 
gives us confidence that our simple model accurately de- 
scribes the system. However, the model is a poorer fit to the 
data after eclipse. This is most likely because the accretion 
disc in SDSS 1507-1-5230 is optically thin, allowing the bright 
spot to remain visible when on the far side of the accretion 
disc. Such an effect is not included in our simple model. To 
limit the effects that the poor fit to post-eclipse data may 
have on the resulting system parameters we excluded the 
regions shown in red (fight grey) in figure [5] from the fit. 

In order to estimate the errors on each parameter once 
the best fit had been found, we perturbed one parameter 
from its best fit value by an arbitrary amount (initially 5 
per cent) and re-optimised the rest of them (holding the 
parameter of interest, and any others originally kept con- 
stant, fixed). We then used a bisection method to deter- 
mine the perturbation necessary to increase by 1, i.e. 

— Xmin — ~ 1- The difference between the perturbed 
and b est-fit values of the parameter gav e the relevant Ict 
error (iLampton. Margon fc Bowveij 1 19761 ) . In the absence 
of systematic errors, or fiaws in the model's assumptions, 
the method described above produc es errors and parame - 
ters which are robust, and unique (|Littlefair et al.ll2006l ). 
We are able to estimate the effects of systematic errors in 
our photometry by fitting each of our eight oberved eclipses 
separately, and comparing the spread of results with the for- 
mal errors produced by the method above. Only the g' - and 
r'-band lightcurves are of sufficient quality to do this. We 
find that the spread in parameters estimated from each indi- 
vidual r'-band eclipse are in good agreement with the formal 
errors, whilst the spread in parameters estimated from the 
individual g'-band eclipses suggest that the formal errors 
underestimate the true uncertainty in some parameters, by 
up to a factor of two. The g'- and r'-band lightcurves have 
almost identical signal/noise ratios, and so we might ex- 
pect them to suffer from systematic errors in the photometry 
at roughly the same level. On the other hand, the g'-band 
lightcurve shows much more flickering prior to eclipse (the 
RMS before eclipse is about 3 times higher in g' than r'). It 
is possible that the presence of flickering affects the model 
estimation in a systematic way. In this case, the formal er- 
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ror is an appropriate estimate. Alternatively, the g'-band 
lightcurve might be more affected by systematic errors in 
the photometry than the r'-band lightcurve, in which case 
the formal errors underestimate the true uncertainty by up 
to a factor of two. Since it is not clear which error estimate 
is more reliable, in each case we adopt the largest of the two 
estimates. These errors are shown in table [S] 

The data were not good enough to determine the white 
dwarf limb-darkening coefficient, f/w, accurately. To find an 
appropriate limb-darkening coefficient, we obtained an es- 
timate of the effective temperature and mass of the white 
dwarf from a first iteration of the method below, and as- 
suming a limb-darkening coefficient of 0.5. The mass and 
effective temperature were t hen used in conjuncti on with 
the stellar atmosphere code of lGansicke et all l| 19951 ) to gen- 
erate angle-dependent white dwarf model spectra. To con- 
vert the spectra to observed fluxes the model spectra were 
folded through passbands corresponding to the instrumen- 
tal response in each fllter; the effects of the SDSS filter 
set, the ULTRACAM CCD responses and the dichroics used 
in the instrumental optics were taken into account. These 
ffuxes were then fit as a function of the limb position in 
order to derive limb-darkening parameters for each band. 
The values adopted are given in table |31 A second itera- 
tion using these values for the limb-darkening parameter 
gave the final values for each parameter. Comparison of 
the limb-darke ning parameters prov i ded by different atmo - 
sphere models (iGansicke et al.lll995l : iHubenv fc LanjflQQSl ) 
suggests an uncertainty in f/w of ~ 5 percent. This corre- 
sponds to an uncertainty in R^/a of ~ 1 percent, which was 
added in quadrature to the error in R^/a provided by the 
fitting process. 

We calculated the remaini ng system parameters foUow- 
ing the method described in iLittlefair et al.l l|2006l ). The 
model fitting provides estimates of q, A</) and Rv,/a. The 
orbital inclination i is determined f rom q and A<^, using 
geometrical arguments l|Bailevl Il979l ). Independent fits to 
the r'-, g'- and w'-band lightcurves yield consistent values 
of q, and R„/a - a weighted mean of these values was 
adopted for the analysis that follows. The model also yields 
fluxes for the white dwarf component, derived by fitting the 
size of the white dwarf ingress and egress. The signal-to- 
noise ratio of our data is good, and the white dwarf fluxes 
are likely to be dominated by systematic errors. We added 
systematic errors of 1 percent to our white dwarf fluxes to 
account for this. A white dwarf temperature was found by 
fltting the white dwarf co lours to the predicti ons of white 
dwarf model atmospheres (|Bergeron et al.|[l995l ) . A mass for 
the white dwarf can then be derived from Kepler's 3'^'' Law, 
the orbital period, the mass ratio and a mass-radius rela- 
tionship for the white dwarf. We adopted a n appropriat e 
mass-radius relationship for the white dwarf l|Woodlll995l ). 
taking into account the effective temperature found above . 
Comparison of different white dwarf m odels l|Woodl [l995l : 
[Bergeron et al. l ll995l : |Panei et al.ll2000l ). differing hydrogen 
envelope masses and correction to different effective temper- 
atures revealed that the dominant source of uncertainty in 
the white dwarf mass is the uncertainty in R„/a. Once the 
white dwarf mass is known, the mass of the donor star fol- 
lows from q. The volume equivalent radius of the donor star 
can be calculated, assuming the donor star fills its Roche 
Lobe. A distance to the system was derived by comparing 
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Orbital phase 

Figure 2. The phased folded u'g'r' (from top to bottom) 
lightcurves of SDSS J1507, fitted separately using the model de- 
scribed in section 13.31 The data (black) are shown with the fit 
(red) overlaid and the residuals plotted below (black). Below are 
the separate lightcurves of the white dwarf (blue), bright spot 
(green), accretion disc (purple) and the secondary star (orange). 
Data points ignored in the fit are shown in red (light grey). 



the white dwarf fluxes in t able |3] and the predicted fluxes 
from [Bergeron et al.l (Il995l ); the uncertainty in the distance 
is dominated by the uncertainty in the white dwarf temper- 
ature. The final adopted system parameters are shown in 
tabled 
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Table 3. Parameters fitted using a modified version of the model of iHorne et alj l|l994l V The fluxes of each component are also shown. 
Prior to fitting, the data were phase-folded and binned by three data points. Note that the orbital inclination i is not a fit parameter 
but is calculated using q and A4>. 



Band u' g' ¥ 

Mass ratio q 0.060 ± 0.002 0.0630 ± 0.001 0.0622 ± 0.001 

Eclipse phase width At/) 0.0405 ± 0.0001 0.04054 ± 0.00005 0.0403 ± 0.0002 

Outer disc radius R^ja 0.343 ± 0.007 0.329 ± 0.004 0.330 ± 0.004 

White dwarf limb darkening C/„ 0.457 ± 0.022 0.346 ± 0.017 0.298 ± 0.015 

White dwarf radius R^ja 0.0177 ± 0.0008 0.0169 ± 0.0006 0.0170 ± 0.0005 

bright-spot scale Sja 0.040 ± 0.008 0.028 ± 0.004 0.032 ± 0.004 

bright-spot orientation 161? ± 2? 157? ± 6? 157? ± 6? 

Isotropic fiux fraction fi^o 0.29 ± 0.02 0.30 ± 0.06 0.31 ± 0.05 

Disc exponent h ± 20 -2.0 ± 2.0 -0.4 ± 0.5 

Phase ofll'set (/)o 38±8xl0-5 15±4xl0-5 10 ± 8 X lO"'^ 

of fit 2555 2108 1470 

Number of datapoints v 1269 1269 1269 

Flux (mjy) 

White dwarf 0.0809 ± 0.0012 0.1260 ± 0.0004 0.1088 ± 0.0004 

Accretion disc 0.0061 ± 0.0022 0.0077 ± 0.0011 0.0182 ± 0.0010 

Secondary 0.0081 ± 0.0013 0.0048 ± 0.0005 0.0058 ± 0.0007 

bright-spot 0.0593 ± 0.0007 0.0651 ± 0.0002 0.0072 ± 0.0003 



Table 4. System param e ters o f SDSS J1507 derived using mass- 
radius relation of IWoodl lll995l ). at the appr opriate Tw. Ry is the 
volume radius of the secondary's Roche lobe llEgeletonlll983l) . The 
weighted means of the appropriate values from Table |3] are used 
for the system parameters. 



Inclination i 
Mass ratio q = M^/M^ 
White dwarf mass M„ /Mq 
Secondary mass M^/Mq, 
White dwarf radius ij„ / Rq 
Secondary radius R^/Rq 
Separation a/RQ 

White dwarf radial velocity K^/hm 
Secondary radial velocity Kr/km 
Outer disc radius R^/a 
White dwarf temperature Tw /K 
Distance d/pc 



83? 63 ±0? 05 
0.0623 ± 0.0007 
0.90 it 0.01 
0.056 ± 0.001 
0.0091 ± 0.0001 
0.096 ±0.001 
0.535 ±0.002 
34.1 ±0.4 
548 ± 2.0 
0.333 ± 0.003 
11000 ±500 
160 ± 10 



4 DISCUSSION 

4.1 The evolutionary status of SDSS J1507 

The orbital period distribution of CVs shows a sharp cut-off 
aro und ~ 80 mins, knowri as the period minimum (see figure 
1 of IPretorius et al1l2007l . for example). The existence of a 
period minimum is easily explained in terms of the response 
of the secondary star to mass loss ljPaczvnskilll98ll ). Above 
the period minimum, the secondary star shrinks in response 
to mass loss and the orbital period decreases. As the sec- 
ondary star nears the substellar limit, however, the thermal 
time-scale becomes longer than the mass-transfer timescale. 
When this happens the secondary is not able to shrink fast 
enough and the period begins to increase. A long-standing 
problem with this theory is that theoretical models persis- 
tently predict a perio d minimum around 67 minutes (e.g. 
iKolb &i Baraffd I1999I ). about 10 minutes short of the ob- 
served period minimum. With an orbital period of 67 min- 



utes and a brown dwarf donor, SDSS J1507 could represent 
the "true" period minimum for CVs. Indeed, the secondary 
star mass and radius presented in table|4]are consistent with 
the predicted mass and r adius for a system nea r the predicted 
mmimum orbital period feolb fc Bara ffelll999l '). There is one 
major problem with this scenario however; where are all the 
systems between 67 minutes and the observed period mini- 
mum at ~ 80 minutes? Population synthesis models predict 
a large nu mber of systems near the period minimum (e.g 
lKolb| ll993^. We can think of no plausible explanation for 
why these systems should be absent and thus conclude that 
the observed period minimum at ~ 78 minutes represents 
the true period minimum for CVs. If this is indeed the case, 
we must then explain why SDSS J1507 has an orbital period 
significantly below 78 minutes. 

Calculations show that CVs containing evolved sec- 
ondary stars, in which a significant fraction of hydrogen in 
the core is processed prior to mass transfer, can reach or- 
bital periods significantly below the period minim um (e.g 
iThorstensen et al.|[20o3 : iPodsiadlowski et al.lliooal ') . Indeed, 
there is good reason to suspect that the two other CVs 
with periods below the observed period minimum contain 
evolved secondary sta rs. EI Psc has a second ary which is 
unusually hot (K4±2 - IThorstensen et al] 2002!) and unusu- 
ally massive (~ O.I2M0 - ISkillman et all [2002') for its or- 
bital period, and the secon dary star in EI Psc has a very 
large N/C abundance ((Gansicke et al.ll2003. 1. The mass and 
spectral type of the secondary star in EI Psc agree well 
with an evolutionary sequence in which a 1.2Mq star begins 
mass transfer close to the end of its hydrogen burning phase 
(|Thorstensen et al]|2002t ) . V485 Cen, too, has a mass which 
is lar ge for its orbital period (i^O.14M0 - lAugusteiin et ahl 
Il996l ). and both stars have He I A6678/Ha ratios higher than 
those typically see n in SU UM a systems (lAugusteiin et ahl 
1996; Thorstcnsen et al]|2002t ). suggestive of an enhanced 
Helium abundance in the donor. 

In the case of SDSS J1507, however, an extremely 
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Figure 3. Model fits to the SDSS optical spectrum of SDSS 
J1507. The light grey curve is the actual spectrum. The white 
dwarf contribution is plotted as a solid red line whose flux in- 
creases towards the blue. The accretion disc contribution is a 
solid green line and the contribution from an M9V star is plotted 
as a solid red line whose flux increases towards the red. 



evolved secondary star cannot be the explanation for the 
short orbital period. The mass of the secondary (Air = 
0.056 ± O.OOIM0) is a factor of two below the secondary 
star masses in EI Psc and V485 Cen. The absence of sec- 
ondary star features in the optical spectrum also suggests 
a rather cool secondary star. We modelled the spectrum 
of SDSS J1507 usi ng three com ponents following the pre- 
scription outlined in lGansicke et al . (2006). The white dwarf 
spectrum is computed for the parameters listed in tabled] 
using SYNSPEC/TLUSTY (|Hubenv fc Lanj|l995l '). The ac- 
cretion disc is des cribed in terms of an isothermal/isobaric 
slab of hydrogen (|Gansicke et al.lll999l ) and the secondary 
star is represented by M /L dwarf spectral templates from 
iBeuermann et al.l (|l998l ) and lKirkpatrick et al.l l|2000l ). The 
secondary star's radius is fixed to 0.097J?© (tabled . Free pa- 
rameters in this model are the temperature, column density 
and flux scaling factor of the hydrogen slab as well as a com- 
mon scaling factor for the white dwarf /secondary star fluxes. 
A plausible flt is shown in flgure|3] which implies a spectral 
type later than M9, which corresponds to an effective tem- 
perature of T^^ ^ < 2450K IVrba et al.ll2004l). This compares 
to a spectral type in EI Psc of K4±2 ( Thorstensen et al] 
|2002| ). Also, the He I A6678/Ha ratio in SDSS J1507 is 0.07. 
This is significantly lower than the ratio of 0.28 measured 
in EI Psc and closer to the average ratio amongst SU UMa 
stars of 0.13 (Thorstensen et al. 2002). The He I A6678/Ha 
ratio depends sensitively on density and temperature as well 
as abundance, so this evidence is not conclusive, but it does 
agree with the general picture that the donor in SDSSJ1507 
is not as evolved as those in EI Psc and V485 Cen. 

There is, of course, a continuum of possibilities between 
the highly evolved secondaries in EI Psc and V485 Cen, and 
a secondary with a normal hydrogen abundance. In order to 
see whether a moderately evolved secondary star could be 
compatible with the observed system parameters we com- 
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Figure 4. The system parameters of SDSS J1507, compared to 
the predictions of some illustrative evolutionary sequences. Solid 
red line: Solar hydrogen abundance, J = 2Jc/j. Blue dashed line: 
an evolved sequence with hydrogen fraction by mass H = 0.47, 
J = 2Jqii. Black dashed line: H = 0.3, J = 2,Jqji. Solid cyan 
line: H = 0.2, J = 2Jgb,- Dashed magenta line: H = 0.3, J = 
Jqr. Red dashed line: sequence with Solar abundance, an initial 
donor mass of O.OTMq and J = 2jQfi where Jqr is the angular 
momentum loss rate predicted by gravitational radiation. 



puted a set of illustrative evo lutionary seque n ces, b ased 
upon the models described in iKolb fc BaraS (|l999l ). To 
ensure our sequences reproduced the observed period mini- 
mum, we adopted an arbitrarily enhanced angular momen- 
tum loss rate, twice that predicted by gravitational radia- 
tion. We also included the effects of distortion of the sec- 
ondary due to t idal a nd rotational forces as quantified by 
iRenvoize et~aLl (|2002h . Models were calculated for an ini- 
tially unevolved secondary star, as well as sequences with 
initial hydrogen fractions (by mass) oi H = 0.47, 0.3 and 
0.2. For H < 0.4, sequences were started with initial donor 
masses below 0.2Mq and with artificially enhanced abun- 
dances. This avoided having to search for the higher-mass 
progenitors of these systems (the progenitors are expected 
to have an intial mass > IMq and to start mass transfer 
toward the end of the central hydrogen burning phase, most 
likely for a central hy drogen msiss frac t ion £ 0.2 according 
to the calculations of iBaraffe fc Kolbl l|2000l ')'). In order to 
quantify the effects of our ad-hoc assumptions about angu- 
lar momentum loss on the evolved sequences, the H — 0.3 
sequence was also calculated with an angular momentum 
loss rate equal to that predicted by gravitational radiation. 
These sequences are shown in figure |4l From this figure it 
is apparent that both the sequence with H = 0.2 and en- 
hanced mass transfer, and the sequence with H = 0.3 and 
"normal" mass transfer are broadly consistent with the ob- 
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served mass, radius and effective temperature limit of SDSS 
J1507. 

There are, however, a number of difficulties encountered 
when we try to explain the observed properties of SDSS 
Jf507 with evolved sequences. The evolved sequences do 
not provide a very good fit to the observed mass and ra- 
dius; if we assume a consistent (enhanced) angular momen- 
tum loss rate for both non-evolved and evolved sequences, 
then the hydrogen fraction must be lower than H <0.2 in 
order to reach orbital periods below ~ 67 minutes. Such 
evolved sequences under-predict the secondary mass. At a 
period of 66 mins, the sequence with H = 0.2 has a mass 
of 0.045M©, which is inconsistent with the observed mass 
of Mr = 0.056 ± O.OOIM©. A good fit can be obtained with 
a sequence with H — 0.3 and angular momentum loss rates 
equal to those predicted by gravitational radiation. Thus, 
to consistently explain both the system parameters of SDSS 
J1507 and the observed period minimum requires that the 
angular momentum loss rate be different in evolved and non- 
evolved sequences. There is no reason why this might not be 
so, but it does imply that the models need some degree of 
"tuning" to consistently explain the secondary star mass in 
SDSS J1507 and the observed period minimum. Addition- 
ally, sequences with an evolved secondary star are only just 
consistent with the upper limit to the effective temperature 
derived above, and fail to explain the slightly lower than 
average He I A6678/H[q ratio in SDSS J1507. It is there- 
fore difficult, though by no means impossible, to explain the 
short period of SDSS J 1507 with an evolved secondary star. 

A second possibility is that SDSS J1507 forme d directly 
from a detached white dwarf/brown dwarf binary. IPolitanol 
showed that CVs forming via this route could have 
orbital periods as short as 46 minutes. iMaxted et al.l (l2006l ) 
recently identified the detached white dwarf /brown dwarf 
binary WDQ137-349, with a white dwarf mass of Q.4M0 and 
a period of 116 minutes, showing that progenitors for such 
CVs can in principle survive the common envelope phase (al- 
though the white dwarf mass is very much lower in WD0137- 
349 than in SDSS J1507). Furthermore, the observed mass 
and radius of the secondary star in SDSS J1507 agree well 
with the predicted values for a system in which the initial 
donor mass was O.O7M0. Given the excellent agreement be- 
tween observed and predicted system parameters, and the 
fact that it is the simplest evolutionary scenario consistent 
with the data, we believe that it is likely that SDSS J1507 
formed directly from a white dwarf/brown dwarf binary. 

4.2 CVs forming with brown dwarf secondaries 

There are now two CVs known which have secondary stars 
with measured masses below the substell ar limit: SDSS 
J1507 and SDSS J103533.03-H055158.4 (|Littlefair et all 
l2006l 'l. With an orbital period of 82.1 minutes, SDSS J1035 is 
consistent with a CV which has evolved past the period min- 
imum. As we argue above however, SDSS J1507 has likely 
formed directly from a detached white dwarf/brown dwarf 

binary. It is the first such system to be ide ntified. 

The puzzle is why there are not more. iPolitancl (|2004l ) 
found that 15% of the zero-age CV population should have 
formed with brown dwarf secondaries and have orbital peri- 
ods below the observed period minimum. Since the evolution 
of these objects is slow, we might expect the present-day CV 



population to be broadly similar in this respect l|Politanol 
2004). Discounting V485 Cen and EI Psc on the basis of 
their evolved secondaries, there is now evidence for at least 
one such system, ou t of the 647 CVs with known periods 
jRitter fc Kolbll2006l) . Could observational selection effects 
explai n the discrepancy? The catalogue of iRitter fc Kolbl 
(|2006l ) is compiled using a wide range of surveys with vary- 
ing selection effects; some CVs are identified using variabil- 
ity, others from their blue colours and others as a result of 
their large outbursts (ICansickdlioOsI ). However, the mass 
transfer rates in CVs forming from brown dwarf secondaries 
are predicted to be comparable to tho se of CVs above the 
period minimum (|Kolb fc Baraffelll 99^1 , of which large num- 
bers exist in the Ritter catalogue. Furthermore, the presence 
of SDSS J1507 amongst the SDSS CV sample suggests that 
this survey, at least, is sensitive to these systems. Within the 
SDSS CV sample there are currently 90 CVs with known pe- 
riods; even within this sub-sample there is clearly a marked 
absence of CVs below the observed period minimum (~ 1% 
compared to 15%). It is therefore unlikely that selection ef- 
fects can explain the absence of significant numbers of CVs 
below the observed period minimum. 



4-2.1 Common envelope efficiency 

One way of reducing the predicted number of CVs forming 
with brown dwarf secondaries is to suggest that they will not 
survive the common envelope phase. If the orbital energy is 
insuffi cient to remove the comm on envelope, a merger will 
occur. iPolitano fc Welled (|2007h showed that, to reconcile 
the observed and predicted number of CVs forming with 
brown dwarf secondaries, the common envelope efficiency 
parameter, Oce, would have to drop to zero for secondary 
masses lower than ~ O.IMq. The existence of SDSS J1507 
argues strongly against this scenario as, by definition, SDSS 
J1507 must have survived its common envelope phase. To 
determine the implications of thi s for common enve l ope e f- 
ficiency, we adopt equation 1 of iPolitano fc Weilerl l|2007h . 
and assume the primary in S DSS J1507 was init ially a main 
sequence star of 3.5-5.5M0 (|Dobbie et al.l [20061 ). We adopt 
a value of 7 = 1 for the dimensionless parameter describing 
the structure of the primary, and assume no mass loss oc- 
curs prior to the common envelope phase. We assume that 
the binary components are co-eval, and determine a lower 
limit to ace by requiring that, at the end of the common 
envelope phase, the brown dwarf remains smaller than its 
Roche Lobe. For a progenitor mass of 3.5M0, we find that 
Qce^ 0.4, whilst for a progenitor mass of 4.5M0, we find 
For larger progenitor masses ace exceeds 1. Sig- 
nificant uncertainties about the evolutionary state of the 
primary at the onset of the common envelope phase means 
that our results are not unambiguous, and so these values 
of Qce should be treated with caution, but our results do 
indicate that SDSS J1507 can survive the common envelope 
phase for physically plausible common envelope efficiencies. 
We note that WD0137-349 provides independent confirma- 
tion that close binaries with brown d warf secondar i es ca n 
survive the common envelope phase ()Maxted et al.ll2006l ). 
It is therefore unlikely that the shortage of CVs forming 
with brown dwarf secondaries is a result of mergers during 
the common envelope phase. 
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^.2.2 The brown dwarf desert 

Whilst the simulations of lPolitano fc Weilej l|2007l ) assume 
a flat mass ratio distribution, observations of companions to 
solar-type stars indicate a relative scarcity of brown dwarf 
companions com pared with either planetary mass or stellar 
companions fe.g. lDuquennov fc Mavor 1991). This has been 
termed the "brown dwarf desert". Icrether fc Lineweaverl 
l|2006t ) found that nearby binaries with solar-type primaries 
and orbital separations ^ 3 AU are ~ 10 times more likely 
to have stellar companions than brown dwarf companions. 
Around 75% of the simulated zero-age CVs forming with 
brown dwarf secondaries had progenitors consisting of solar- 
type primaries with orb ital separations less than 3 AU 
(Politano fc Weiler ''2007''). The dearth of CVs with periods 
below the period minimum could therefore be a simple con- 
sequence of a shortage of progenitors for these CVs. Since 
the observed properties of SDSS J 1507 seem to indicate that 
the shortage of CVs with periods below the period minimum 
cannot be explained by observational selection effects or by 
mergers during the common envelope phase, we conclude 
that the dearth of CVs with periods shorter than the pe- 
riod minimum provides independent support for the reality 
of the brown dwarf desert. 



limit, making it the second confirmed substellar donor in a 
cataclysmic variable. 

We attempt to explain the short orbital period of SDSS 
J1507 as a result of nuclear evolution of the secondary star. 
We argue that given the secondary star's low mass, the lack 
of any visible signature of the secondary star and the low 
He I A6678/HQf line ratio, together with the difficulty in con- 
sistently explaining the secondary star mass in SDSS J1507 
and the observed minimum period, argue against this sce- 
nario. Instead, we suggest that SDSS J1507 may have formed 
directly from a detached white dwarf/brown dwarf binary. 
If this scenario is correct, SDSS J1507 will be the first CV 
known to have formed this way. The effective temperature of 
the secondary, as well as its surface abundances, ought to be 
quite different in the two cases; deep infrared spectroscopy 
of SDSS J1507 would be highly desirable. 

If SDSS J1507 did indeed form directly from a detached 
white dwarf/brown dwarf binary, its observed properties ar- 
gue against the scarcity of such systems being explained by 
observational selection effects, or mergers during the com- 
mon envelope phase. The best remaining explanation is that 
the progenitors of CVs are themselves scarce, providing in- 
dependent evidence for the existence of the brown dwarf 
desert. 



4.3 The white dwarf temperature and ZZ Ceti 
pulsations 

The observational ZZ Ceti instability stri p ranges in effec- 
tive temperature from 12270 to 10850K (|Gianninas et al.] 
I2OO5I ). Current theoretical models predict that the ZZ Ceti 
instability strip should be pure: all the isolated DA white 
dwarf stars, with temperatures within the instability strip 
limits should pul sate, with a small dependency on mass 
l| Arras et aLlbOOeD . With an effective temperature of llOOOi 
500K, the white dwarf in SDSS J1507 is inside the conven- 
tional ZZ Ceti strip. Lomb-Scargle periodograms reveal no 
evidence for ZZ Ceti pulsations, however, despite our rapid 
photometry and good signal-noise ratio. 

In fact, almost all of the observed pulsating white 
dwarfs in CVs are outsid e the conventional ZZ Ceti s trip 
l|Szkodv et al.l |2007| . |2002| : lArauio-Betancor et al.l 120051 ') . a 
result which is probably explained by the effects of a spread 
in white dwarf masses and accre tion of small amo unts of 
Helium from the secondary star (| Arras et al.ll200^ . These 
effects may explain the absence of ZZ Ceti pulsations in 
SDSS J1507. 



5 CONCLUSIONS 

We present high-speed, three-colour photometry of the 
eclipsing cataclysmic variable SDSS J150722.30-I-523039.8. 
We measure an orbital period of 66.61 minutes, placing 
it below the observed "period minimum" for cataclysmic 
variables. By fitting a parameterised model to the eclipse 
lightcurves we obtain a mass ratio of 5 = 0.0623 ± 0.0007, a 
primary mass of Mw = 0.90 ± O.OIMq and a secondary 
mass and radius of Mr = 0.056 ± O.OOIMq and Ry = 
0.096 ± O.OOli?0, respectively. The secondary star in SDSS 
J1507 has a mass substantially below the hydrogen burning 
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